Dynamics and length distribution of microtubules under force and confinement 
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We investigate the microtubule polymerization dynamics with catastrophe and rescue events for 
three different confinement scenarios, which mimic typical cellular environments: (i) The micro- 
tubule is confined by rigid and fixed walls, (ii) it grows under constant force, and (iii) it grows 
against an elastic obstacle with a linearly increasing force. We use realistic catastrophe models 
and analyze the microtubule dynamics, the resulting microtubule length distributions, and force 
generation by stochastic and mean field calculations; in addition, we perform stochastic simulations. 
Freely growing microtubules exhibit a phase of bounded growth with finite microtubule length and 
a phase of unbounded growth. The main results for the three confinement scenarios are as follows: 
(i) In confinement by fixed rigid walls, we find exponentially decreasing or increasing stationary mi- 
crotubule length distributions instead of bounded or unbounded phases, respectively. We introduce 
a realistic model for wall-induced catastrophes and investigate the behavior of the average length 
as a function of microtubule growth parameters, (ii) Under a constant force the boundary between 
bounded and unbounded growth is shifted to higher tubulin concentrations and rescue rates. The 
critical force f c for the transition from unbounded to bounded growth increases logarithmically with 
tubulin concentration and the rescue rate, and it is smaller than the stall force, (iii) For microtubule 
growth against an elastic obstacle, the microtubule length and polymerization force can be regulated 
by microtubule growth parameters. For zero rescue rate, we find that the average polymerization 
force depends logarithmically on the tubulin concentration and is always smaller than the stall force 
in the absence of catastrophes and rescues. For a non-zero rescue rate, we find a sharply peaked 
steady-state length distribution, which is tightly controlled by microtubule growth parameters. The 
corresponding average microtubule length self-organizes such that the average polymerization force 
equals the critical force f c for the transition from unbounded to bounded growth. We also investigate 
the force dynamics if growth parameters are perturbed in dilution experiments. Finally, we show 
the robustness of our results against changes of catastrophe models and load distribution factors. 

PACS numbers: 87.16.Ka, 87.16.-b 
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I. INTRODUCTION 



Microtubules (MTs) are one of the main components 
of the cytoskeleton in eukaryotic cells. Their static and 
dynamic properties are essential for many cellular pro- 
cesses. MTs serve as pathways for molecular motor pro- 
teins [1] and contribute to cell stiffness [2] . Dynamic MTs 
play a crucial role in the constant reorganization of the 
cytoskeleton, and single MTs can generate polymeriza- 
tion forces up to several pN [3]. These forces are used 
in various intracellular positioning processes [4], such as 
positioning of the cell nucleus [5] or chromosomes during 
mitosis , establishing cell polarity [6] , or regulation of cell 
shapes [7, 8]. In many cellular processes, MTs establish 
and maintain a characteristic length in response to forces 
exerted, for example, from the confining cell cortex [7]. 

The fast spatial reorganization of MTs is based on the 
dynamic instability: Polymerization phases are stochasti- 
cally interrupted by catastrophes which initiate phases of 
fast depolymerization; fast depolymerization terminates 
stochastically in a rescue event followed again by a poly- 
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merization phase [9]. This complex dynamic behavior 
with catastrophes and rescue events is central to a rapid 
remodelling of MTs in the cytoskeleton, but it also affects 
their ability to generate polymerization forces. We will 
show that, in general, the dynamic instability decreases 
the average polymerization force of a single MT. 

In this article we theoretically investigate the polymer- 
ization dynamics of a single MT under force or confine- 
ment and in the presence of the MT dynamic instability. 
We use a coarse-grained polymerization model with dy- 
namic instability and characterize spatial and temporal 
behavior in three different scenarios, which mimic typ- 
ical cellular environments that can also be reproduced 
in vitro: (i) Confinement: The MT is confined between 
fixed rigid walls, which cannot be deformed by the micro- 
tubule. Such confinement is realized in fixed solid cham- 
bers [10]. (ii) Constant force: A constant force is acting 
on the MT. Constant forces can be realized by optical 
tweezers with a force clamp control [11]. (iii) Elastic ob- 
stacle: The microtubule grows against an elastic obstacle, 
which resists further growth by a force growing linearly 
with displacement. Elastic forces can be realized by op- 
tical tweezers without force clamp [12, 13]. For all three 
confinement scenarios (i)-(iii) we focus on the resulting 
length distributions of MTs and for scenarios (ii) and 
(iii), we calculate the polymerization force that a single 
MT can generate. 
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Dynamic MTs also initiate regulation processes or are 
subject to regulation. Dynamic MTs can activate or de- 
activate proteins upon contacting the cell membrane [14] , 
or they can activate actin polymerization within the cell 
cortex [15, 16]. At the same time, polymerizing MTs are 
also targets of cellular regulation mechanisms [17], which 
affect their dynamic properties. 

The dynamic instability of MTs enables various regu- 
lation mechanisms of MT dynamics. Catastrophes and 
rescues result from the hydrolysis of GTP-tubulin within 
MTs. When GTP-tubulin is incorporated into the tip 
of a growing MT, it forms a stabilizing GTP-cap. The 
loss of this GTP-cap due to hydrolysis of GTP-tubulin 
to GDP-tubulin, causes a catastrophe [9]. In living 
cells, there are various microtubule associated proteins 
(MAPs) that either stabilize or destabilize microtubules 
and regulate microtubule dynamics both spatially and 
temporally [18]. Recently, the importance of MAPs 
in connection with the plus end of growing MTs has 
been recognized [19]. Stabilizing MAPs bind to assem- 
bled MTs, thereby reducing the catastrophe rate or in- 
creasing the rescue rate. Destabilizing MAPs such as 
OP18/stathmin bind to GTP-tubulin dimers, thus de- 
creasing the available GTP-tubulin concentration, which 
in turn decreases the growth velocity of the GTP-cap and 
makes catastrophes more likely. Therefore, such mecha- 
nisms can regulate basic parameters in our model, such 
as the available GTP-tubulin concentration or the rescue 
rate, and we will systematically study their influence on 
the generated polymerization force for the three confine- 
ment scenarios (i)-(iii). 

The paper is structured as follows: In Sec. II, the MT 
model and the basic notation are introduced. We also 
discuss the catastrophe model and the underlying hy- 
drolysis mechanism in the absence and in the presence 
of a resisting force. Section III deals with the simulation 
model. In Sees. IV, V, and VI, results for the three sce- 
narios (i)-(iii) are presented and discussed. In Sec. VII, 
the elastic obstacle is reconsidered using an alternative 
catastrophe model based on experimental measurements 
to show that our results are robust with respect to this 
change in the catastrophe model. In Sec. VIII, we show 
that our results are also robust with respect to possible 
generalization of the force-velocity relation by introduc- 
ing load-distribution factors. Section IX contains a final 
discussion and outlook. 



II. MICROTUBULE MODEL 

A. Single MT dynamics 

The MT dynamics in the presence of its dynamic in- 
stability is described in terms of probability densities and 
switching rates [20, 21]. In the growing state, a MT poly- 
merizes with average velocity v + . The MT stochastically 
switches from a state of growth (+) to a state of shrink- 
age (— ) with the catastrophe rate w c . In the shrinking 



state, it rapidly dcpolymerizes with an average velocity 
V- ~ 3 • 10 _7 m/s (Table I). With the rescue rate u> r 
the MT stochastically switches from a state of shrinkage 
back to a state of growth. We model catastrophes and 
rescues as Poisson processes such that (t+) = 1/uj c and 
(t_) = l/uj r are the average times spent in the growing 
and shrinking states, respectively. The stochastic time 
evolution of an ensemble of independent MTs, growing 
along the x-axis, can be described by two coupled mas- 
ter equations for the probabilities p+(x, t) and p-(x,t) 
of finding a MT with length x at time t in a growing or 
shrinking state, 

d t p+{x,t) = -ui c p + (x,t) + u r p_(x,t) -v + d x p + (x,t) 



d t p-(x,t) = (J c p + (x,t) - u r p-(x,t) + v-d x p-(x,t). 



(1) 
'(2) 



In the following, we will always use a reflecting boundary 
at x = 0: A MT shrinking back to zero length undergoes 
a forced rescue instantaneously. This corresponds to 



v+p+(0,t)=v-p-(0,t). 



(3) 



A more refined model including a nucleating state has 
been considered in Ref. [22]. For a constant and fixed 
catastrophe rate uj c , eqs. (1) and (2) together with the 
boundary condition (3) can be solved analytically on the 
half-space x > 0, and we can determine the overall prob- 
ability density function (OPDF) of finding a MT with 
length x at time t, P(x,t) =p+(x,t) +p-(x,t) [20, 21]. 
The solution exhibits two different growth phases: a 
phase of bounded growth and a phase of unbounded 
growth. 

In the phase of bounded growth the average length 
loss during a period of shrinkage, u_(r_) = V-/uj r , ex- 
ceeds the average length gain during a period of growth, 
v + (t + ) = v + /u> c . The steady-state solution of P{x,t) 
assumes a simple exponential form P(x) = | A | 1 e ' A ' 
with an average length (x) = |A| and a characteristic 
length parameter 



A 



V+V- 



V+LJ r — V^UJ C 



(4) 



with A -1 < for bounded growth [20]. The transition 
to the regime of unbounded growth takes place at A -1 = 
0, where the average length gain during growth equals 
exactly the average length loss during shrinkage, 



v+ui r 



(5) 



such that (x) diverges. 

In the regime of unbounded growth (A > 0), the av- 
erage length gain during a period of growth is larger 
than the average length loss during a period of shrink- 
age. There is no steady state solution, and for long times 
P(x,t) asymptotically approaches a Gaussian distribu- 
tion [20] 



P(x,t) 



1 



2Vt^DjI 



exp 



(a; - J if 
4Djt 



(6) 
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centered on an average length which approaches linear 
growth (x) ps Jt with a mean velocity J and with diffu- 



sively growing width (x 



2Djt with a diffusion 



constant Dj. The average growth velocity is given by 



V-LO c 



(7) 



because the asymptotic probabilities to be in a growing or 
shrinking state are 7r + = uj r /(uj c +uj r ) and 7r_ = w c /(w c + 
u) r ), respectively. The diffusion constant Dj is 



Dj = 



LO c L0 r (V-{ 



{uj c + uj r y 



(8) 



The transition between the two growth phases can be 
achieved by changing one of the four parameters of MT 
growth, u> c , uj r , or u_. In the following, we will 
use catastrophe models, where the catastrophe rate w c 
is a function of the growth velocity v+, which in turn 
is determined by the GTP-tubulin concentration via the 
GTP-tubulin on-rate w n (assuming a fixed off-rate 0j o f[). 
Moreover, experimental data suggest that V- is fixed to 
values close to ~ 10~ 7 m/s (Table I). As a consequence, 
there are two tunable control parameters left, the GTP- 
tubulin concentration or, equivalently, the tubulin on- 
rate Won and the rescue rate u r . These are the control 
parameters we will explore for MTs in confinements and 
under force. These parameters are also targets for regu- 
lation by MAPs, such as OP18/stathmin, which reduces 
Won by binding to GTP-tubulin dimers or MAP4, which 
increases the rescue rate ui r - 



B. Force-dependent catastrophe rate 

In a growing state, GTP-tubulin dimers are attached 
to any of the 13 protofilaments with the rate w on , which 
is directly related to the GTP concentration. We explore 
a regime u) on = 30, 100s _1 , see Table I. GTP-tubulin 
dimers are detached with the rate w g' = 6s _1 [23] such 
that we can typically assume uj oa 3> oj Q g. In the absence 
of force or restricting boundaries, the velocity of growth 
is given by 



V + (0) = d (w on - M 



(9) 



Here d denotes the effective dimer size d w 8nm/13 ps 
0.6 nm. 

The classical view of the MT catastrophe mechanism is 
based on a purely chemical model of catastrophes, where 
the catastrophe rate u c is determined by the hydroly- 
sis dynamics of GTP-tubulin [9]. When GTP-tubulin 
is incorporated into the tip of a growing MT, it forms 
a stabilizing GTP-cap. In a chemical model, the loss of 
this GTP-cap due to hydrolysis of GTP-tubulin to GDP- 
tubulin directly causes a catastrophe. However, recent 
research indicates that the "structural plasticity" of the 
MT lattice can play a role for the kinetics of catastrophes 
[24] . This structural plasticity mechanism is based on the 



assumption that GDP-tubulin prefers a curved configu- 
ration, which generates additional mechanical stresses in 
the MT by hydrolysis. Also in the presence of structural 
plasticity, the loss of the GTP-cap has a destabilizing 
effect, but the kinetics leading to a catastrophe can be 
more complicated because the initiation of a catastro- 
phe event is similar to the nucleation of a crack in the 
stressed MT lattice within this model. In this article, we 
focus on purely chemical catastrophe models and neglect 
mechanical effects on the catastrophe kinetics. 

Within a chemical catastrophe model, the loss of the 
GTP-cap due to hydrolysis of GTP-tubulin to GDP- 
tubulin triggers a catastrophe immediately. Therefore, 
the catastrophe rate u> c is given by the first-passage 
rate to a state with vanishing GTP-cap and has been 
discussed within a model with cooperative hydrolysis 
[25, 26], where GTP-tubulin is hydrolyzed by a com- 
bination of random and vectorial mechanisms; similar 
models have also been discussed for hydrolysis in F- 
actin [27, 28]. In random hydrolysis, GTP-tubulin is 
hydrolyzed at a random site within the GTP-cap with 
a rate per length r ±2 3.7 • 10 6 m _1 s _1 , while in vecto- 
rial hydrolysis, only GTP-tubulin with adjacent GDP- 
tubulin is hydrolyzed. This results in hydrolysis fronts 
propagating through the microtubule with average ve- 
locity Vh — 4.2 • 10~ 9 m/s. The inverse catastrophe rate 
can then be calculated as the mean first-passage time to a 
state with zero cap length, as a function of hydrolysis pa- 
rameters and v+. With v = v+ — Vh, D = 0.5d(v + + Vh) 
and 7 = 0.5u_D 1 / 3 r~ 1 / 3 the exact analytical result for 
the dimcnsionless catastrophe rate a = Lu c D~ 1 / z r~ 2 /' i is 
given by the smallest solution of 



Ai' ('y 2 — a) = —jAi(~/ 2 — a). 



(10) 



Here Ai denotes the first Airy function and Ai' its deriva- 
tive [29]. We solved eq. (10) numerically and obtained a 
high order polynomial for the function a = a (7). This 
polynomial is used in simulations and analytical calcula- 
tions to compute the catastrophe rate oj c = a(7)£> 1 / 3 r 2 / 3 
as a function of the growth velocity v + , while the hydrol- 
ysis parameters Vh and r are fixed. 

Under a force F, the tubulin on-rate w on is modified 
by an additional Boltzmann factor [30] and the force- 
dependent growth velocity becomes 



v+(F) = d [won exp (-Fd/k B T) 



Woff 



(11) 



Here Fd is the work that has to be done against the 
force F to incorporate a single dimer of size d; ks is the 
Boltzmann constant and T = 300 K the temperature. In 
the following we use the dimcnsionless force 



f = F/F with F = k B T/d, 



112) 



in terms of which the force-dependent growth velocity is 
given by 



v+(f) = d[ 



w on e / - Woff] 



(13) 
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The characteristic force Fq has a value Fq = ksT ' Jd « 
7pN. The dimensionless stall force 

/stall = In (cj on /u; off ) (14) 

is defined by the condition of vanishing growth velocity 
t>+ (/stall) = 0. We typically have / BtaU ~ 1.5,..., 3 or 
F staU 10, ...,20pN for w on = 30, 100 s" 1 . The stall 
force is the maximal force that the MT can generate in 
the absence of catastrophes. We will investigate how the 
forces that can be generated in the presence of catastro- 
phes compare to this stall force. 

The velocity-dependence of the catastrophe rate as cal- 
culated from eq. (10) gives rise to a force-dependence 
lo c = lo c (d_|_(/)). Wc assume that this is the only effect 
of force on the catastrophe rate [39]. As a result, the 
catastrophe rate increases exponentially, when «+(/) is 
decreased by applying a force /, but a finite value is main- 
tained at v + (f) = 0, which is uj c (v + = 0) « 2.9s -1 . Wc 
assume that t>_ is independent of force. For qualitative 
approximations, the force-dependence of the catastrophe 
rate can be described by an exponential increase above 
the characteristic force Fq, 



(15) 



In Sec. VII we introduce an alternative catastrophe model 
which is based on experimental measurements. The ex- 
ponential approximation (15) applies to the catastrophe 
model described above as well as to the alternative catas- 
trophe model, see Fig. 11. Our results are robust for all 
catastrophe models with an exponential increase above 
the characteristic force Fq. Our results do not directly 
apply to more elaborate multi-step catastrophe models 
with more than two MT states [31]. 



III. SIMULATION MODEL 

In the simulations we solve the stochastic Langevin- 
like equations of motion for the length x(t) of a single 
MT using numerical integration with fixed time steps At 
and including stochastic switching between growth and 
shrinkage. In a growing state x(t) is increased by v+At, 
while in a state of shrinkage, it is decreased by V-At. 
In the growing state, v + is calculated from eq. (9) for 
zero force and from eq. (11) under force. In each time 
step a uniformly distributed random number £ e [0, 1] 
is compared to uj r ^ c At. If £ < oj r ^ c At the MT changes 
its state of growth. The catastrophe rate oj c is calculated 
from the high order polynomial obtained from eq. (10) as 
mentioned above. To assure Lu r ^ c At < 1 we used a time 
step At = 0.1s. During the simulations all parameters 
of growth, d = 0.6 nm, r = 3.6 • 10 6 m -1 s _1 , Vh = 4.2 • 
10~ 9 m/s, k B = 1.38 • 10- 23 J/K, T = 300 K, and w off = 
6s , are fixed, see Table II, except for cj on , which is 
varied in the range u> on = 30 — 100s _1 , and u r) which 
is varied in a range u> r = 0.03 — 0.2s -1 , see Table I. 
Averages are taken over many realizations of stochastic 
trajectories. 
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FIG. 1: (a): Schematic representation of the confinement 
and possible MT configurations. From top to bottom: A MT 
growing with v + ; MT shrinks with V-. MT in a state of 
growth and stuck to the boundary wall with i>+ = and lo c ,l- 
(b): Schematic representation of a single MT growing against 
the elastic obstacle. From top to bottom: MT shrinks with 
V-. MT under force F(x) = k(x — x ) with f(x) = F(x)/F , 
v+[f(x)], and force-dependent catastrophe rate u c [f(x)]. 



IV. CONFINEMENT BETWEEN FIXED RIGID 
WALLS 

A single MT is confined to a one-dimensional box of 
fixed length L with rigid boundary walls at x = and 
x = L as shown schematically in Fig. 1(a) [32, 33]. There 
is no force acting on the MT but within the box catastro- 
phes are induced upon hitting the rigid walls. We propose 
the following mechanism for these wall-induced catastro- 
phes: When the MT hits the boundary at x = L, its 
growth velocity v + has to reduce to zero, which leads to 
an increase of the catastrophe rate to uj c x = w c (v+ = 0). 
Since u} c ,l is finite, wall-induced catastrophes are not 
instantaneous but the MT sticks for an average time 
l/w Ci i to the boundary before the catastrophe, which 
is in contrast to previous studies [34]. For the average 
time spent at the boundary before a catastrophe, we find 
wi w 0.29s. The catastrophe rate at the wall, <jJ c ,l, is 
much higher than the bulk catastrophe rate u) c (v+). For 
Won = 50s _1 wc find Uc.l/lOc ~ 2300. 

To include the mechanism of wall- induced catastrophes 
into the description by master equations, we introduce 
the probabilities Q+ and Q_ of finding the MT stuck to 
the boundary in a growing state and in a shrinking state, 
respectively. The stochastic time evolution of Q+ (i) and 
Q-(t) is given by: 

d t Q+{t) = -Lo c . L Q+(t) + u r Q-(t) +v+p+(L) (16) 

d t Q-(t) = +u c , L Q+{t) - w r Q-(*) - \Q-{t). (17) 

The quantity v+p.\-{L) is the flow of probability from the 
interior of the confining box onto its boundary and is 
given by the solution of eq. (1) and (2) for x = L, while 
(v- 1 A)Q_ is the probability current from the boundary 
back into the interior, where A denotes a small interval in 
which the flow i>_Q_ can be measured. This implies that 
there is a boundary condition V-P—(L,t) = {v-/A)Q- 
for the backward current density at x = L, in addition 



5 



to the reflecting boundary condition (3) at x = 0. An 
identical model for wall-induced catastrophes has been 
introduced in Ref. [22] recently. 

In the steady state and in the limit A « we find 

v + 



o, 



-P+(L) 



(18) 
(19) 



and V-P-(L,t) = (u_/A)Q_ = v+p+(L). Eq. (18) 
shows that there is a non-zero probability Q + of find- 
ing a MT in a state of growth and stuck to the bound- 
ary, which is given by the flow of probability from the 
interior of the confining box onto its boundary divided 
by the average time being stuck to the boundary. In 
contrast, eq. (19) states that there is no MT in a shrink- 
ing state and stuck to the wall. This is intuitively clear 
since a MT undergoing a catastrophe begins to shrink 
instantaneously. In the steady state, we solve eqs. (1), 
(2) and (18) simultaneously with the additional nor- 
malization Jq(p + (x) + p-(x))dx + Q + = 1. We find 
v+p+{x) = V-p~(x) and 



P{x) = Ne x/X ^1 + ^ 



(20) 
(21) 



AT 



with A from eq. (4) and a normalization 

-i = \(l + ^)(e L / x -l)+^e L /\ (22) 

\ V-J \ J UJ C ,L 

Equation (20) shows that we find an exponential OPDF 
P{x) in confinement with the same characteristic length 
|A|. If the growth is unbounded in the absence of con- 
finement, which corresponds to A -1 > 0, the OPDF is 
exponentially increasing; if the growth is bounded in the 
absence of confinement, which corresponds to A -1 < 0, 
the OPDF remains exponentially decreasing in confine- 
ment. The same result has been obtained in Ref. [34] 
within a discrete growth model. In independent in vivo 
experiments, both exponentially increasing [35] and ex- 
ponentially decreasing OPDFs [20] have been found. 

In the following we focus on the case A^ 1 > of ex- 
ponentially increasing OPDFs. In the steady state, the 
average length of a MT within the confining box is given 
by 



/ xP(x)dx + Q + L 
Jo 



N i ( 1 + — ) A 2 



(23) 



In the limit of instantaneous wall-induced catastrophes, 
Q + si 0, we obtain 



1 



A 



l-e- L /A V 



i.e., the average MT length (x)/L depends on the two 
control parameters u r and cj on only via the ratio L/X. 
This scaling property is lost if wall-induced catastrophes 
are not instantaneous because eq. (23) then exhibits addi- 
tional v+- and thus w on -dependencies. Within our model 
the increased catastrophe rate at the boundary gives rise 
to an increased overall average catastrophe rate 



Wc.off = w c (v+) + Q + (u c x - lo c (v + )), 



(25) 



(24) 



for which we find w c ,eff — 0.03 s _1 for L = 1 /im and 
w Cje ff ~ 0.006 s -1 for L = 10 /iin as compared to lo c ~ 
0.0015 s -1 for these conditions. 

We set the length of the confining box to L = 1 /xm 
and L = 10 /im, which are typical length scales in ex- 
periments [10, 11] and cellular environments [5], and we 
calculate (x) and Q+ as functions of oj on and ui r . The 
parameter regimes displayed in Figs. 2 and 3 correspond 
to regimes L/X ^> 1 for L = 10/im and L/X <C 1 for 
L = 1 ^m. Results obtained from stochastic simulations 
agree with analytical findings (Figs. 2 and 3). It is clearly 
visible that the size L of the confinement has a significant 
influence on (x), mainly via the ratio L/X. 

The probability Q + to find the MT at the wall in- 
creases with increasing rates in the range of Q + w 
0, ...,0.03 and exhibits only a weak dependency on L, 
see Figs. 3. Even for maximum rates, the probability of 
finding a MT in a growing state and stuck to the wall 
is limited to several percent, due to the large catastro- 
phe rate lu c ,l at x = L. Therefore, in most cases wall- 
induced catastrophes can be viewed as instantaneous, 
and the approximation (24) works well. For increasing 
on-rate w on or rescue rate oj r , the ratio L/X approaches 
L/X ~ Lu r /v- from below. According to the approx- 
imation (24), the mean length (x) then increases and 
approaches (x)/L rj 1/(1 — e~ LUr ^ v ~) — V-/Lui T from 
below. For L = 10 fim, we have t/A > 1 and the 
length distribution is exponential, P(x) ~ e 2 / A . The 
ratio (x)/L saturates at a high value (x) / L w 0.7, ...,0.9 
(Figs. 2 (a),(c)). For L/X > 1 the MT length distribu- 
tion becomes very narrow around the maximal length L. 
In contrast, for L = 1 /im, we have L/X <C 1, and L is 
too small to establish the characteristic exponential de- 
cay of the length distribution. The length distribution 
P{x) is almost uniform, and the ratio (x) /L m 0.5, 0.6 
deviates only slightly from the result (x)/L = 1/2 charac- 
teristic for a broad uniform distribution (Figs. 2(b), (d)). 



V. CONSTANT FORCE 

In the second scenario a constant force F is applied 
to the MT and the right boundary is removed, so that 
the MT is allowed to grow on x £ [0,oo[. According 
to eq. (13) the growth velocity under force is smaller, 
but it remains constant for fixed /. With eq. (10) this 
results in a higher, but also constant, catastrophe rate 
u> c [v+(f)] > u! c [v+(0)]. Since u_ and uj r arc independent 
of force, the stochastic dynamics of the MT is described 
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FIG. 2: The average length (x) as a function of uion and uj r for 
confinement by fixed rigid walls. Data points are results from 
stochastic simulations, lines are analytical results (23). Top 
row: The average length (x) as a function of LJ on for different 
values of u) r = 0.03 s -1 (E]), 0.05 s _1 (B), 0.1 S _1 (0), 0.2 s" 1 (•) 
and 0.3 s~ (A), (a) L = 10 fjxa. (b) L = 1/im. Lower row: 
The average length {x) as a function of uj r for different values 
of w on = 25s- 1 (B),50s- 1 (Q),75s- 1 (A),100s- 1 (E]). (c) L = 
10 /im. (d) L = 1 fixa. 
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FIG. 3: The probability Q+ to find the MT at the 
wall as a function of Lu on and ui r for confinement by 
fixed rigid walls. Data points are results from stochas- 
tic simulations, lines are analytical results (18). Top 
row: Q+ as a function of u) on for different values 
of uj r = O.O3s- 1 (Q),O.O5s- 1 (B),O.ls- 1 (0),O.2s- 1 (») and 
0.3s" 1 (A), (a) L = 10 /im. (b) L = 1 /xm. 
Lower row: Q+ as a function of uj r for different values of 
cjon = 25s- 1 (B),50s- 1 (Q),75s~ 1 (A),100s~ 1 (Q). (c) L = 
10 /im. (d) L = 1 fim. 



by cq. (1) and (2) with the same solutions P(x,t) as in 
the absence of force, but with a decreased velocity of 
growth t>+(/) and an increased catastrophe rate u> c (f) 
[20, 21]. In particular, we still find two regimes, a regime 
of bounded growth and a regime of unbounded growth. 

In the regime of bounded growth P(x,t) is again ex- 
ponentially decreasing, and the force-dependent average 
length is (x(f)) = |A(/)| with the corresponding force- 
dependent length parameter 

v+(f)v- 



A(/) 



(26) 



as compared to eq. (4) in the absence of force. In 
the regime of unbounded growth (x(f)) increases lin- 
early in time with the force-dependent mean velocity 
J(f) = [v + (f)oj r -V-Lu c (f)}/[uj r +GJ c (f)}, cf. eq. (7). The 
MT length distribution P(x, t) assumes again a Gaussian 
form (6) where also the diffusion constant Dj(f) follows 
the same eq. (8) with force-dependent growth velocity 
v + (f) and catastrophe rate w c (/). 

In the presence of a constant force /, the transition 
between bounded and unbounded growth is governed by 
the force-dependent parameter A(/). The regimes of 
bounded and unbounded growth are now separated by 
the condition A~ 1 (f) = 0, which is shifted compared to 
the case / = 0, see Fig. 4(a). The inverse length parame- 
ter A -1 (/) is a monotonously decreasing function of force 
/ and changes sign from positive to negative values for 
increasing force /. Therefore A _1 (/ c ) = or 



v+(fc)u r = u_w c (/ c ), 



(27) 



defines a critical force for the transition from unbounded 
to bounded growth. A single MT exhibiting unbounded 
growth (A _1 (0) > 0) in the absence of force undergoes 
a transition to bounded growth with A _1 (/) < by ap- 
plying a supercritical force / > f c . On the other hand, 
starting with a combination of on-rate w on and rescue 
rate to r and a force /, which results in bounded growth 
with A~ 1 (f) < 0, the MT can still enter the regime of 
unbounded growth by increasing cj on or uj r so that the 
force / becomes subcritical, A _1 (/) > or / < f c . 

Rewriting condition (27) as u + (/ c ) = v-uj c (f c )/uj r > 
and using that v+(f) decreases with /, it follows that the 
critical force is always smaller than the stall force, f c < 
/stall, which satisfies i>+(/ s taii) = 0, and it approaches 
the stall force only for vanishing catastrophe rate. Qual- 
itatively, we can obtain an explicit result for the critical 
force f c by using the approximations of an exponentially 
decreasing growth velocity, v+(f) ~ v + (0)e~f, which is 
valid for uj on ^> u; ff (see eq. (13)), and an exponentially 
increasing catastrophe rate above the characteristic force 
Fq, eq. (15), in the condition (27) for the critical force. 
This leads to 



fc 



2 



v+ (0)w r 



Im 

2 



,diO r 



(28) 



which shows that the critical force grows approximately 
logarithmically with on-rate uj on (note that the catastro- 
phe rate in the absence of force decreases with uj on as 
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FIG. 4: (a): Phase boundary between bounded 

(B) and unbounded growth (UB) as a function of ui on 
and uj r for MT growth under constant force. Data 
points for / = 0(H), 0.3(B), 0.6(0), 0.9(»), 1.2(A), 1.4(A) 
represent results from simulations, lines represent solu- 
tions of v+(f)td r = v -^c(f) for a constant force /. 
(b): Critical force f c as a function of uj on for uj r — 
O.O3s~ 1 (H),O.O5s- 1 (B),O.ls- 1 (0),O.2s~ 1 (»). Data points 
represent results from simulations, lines represent the solu- 
tion of eq. (27) for a fixed combination of Lj on and uj r . 



w c (0) oc l/u> on [26]) and rescue rate u r . A negative f c 
for small on-rates and rescue rates signals that the MT is 
for all forces / > in the bound phase. In Fig. 4(b) we 
show exact results for the critical force f c as a function 
of the on-rate uj on and for different rescue rates ui r from 
solving condition (28) numerically and from stochastic 
simulations. Agreement between both methods is good. 

The condition \~ 1 (f) — specifies the boundary be- 
tween bounded and unbounded growth at a given force /. 
In Fig. 4(a), the resulting phase boundary is shown as a 
function of w on and uj r . There is good agreement between 
numerical solutions of A~ 1 (f) = and stochastic simula- 
tions. With increasing force, the boundary between the 
two regimes of growth shifts to higher values of oj on and 
u> r , and forces up to F ~ 1.4 • Fq can be overcome by a 
single MT in the parameter regimes of w on and uj r con- 
sidered. 



VI. ELASTIC FORCE 

In the third scenario, an elastically coupled barrier is 
placed in front of the MT as shown in Fig. 1(b), which 
models the optical traps used in Rcfs. [11, 13] or the elas- 
tic cell cortex in vivo. If the barrier is displaced from its 
equilibrium position xq by the growing MT with length 
x > .To, it causes a force F(x) = k(x—xo) resisting further 
growth. For x < xq there is no force. We use xq = /jm 
in the case of vanishing rescue rate and xq = 10 /jm in 
the case of finite rescue rate and a spring constant k in 
the range 10" 7 N/m (soft) to lCT 5 N/m (stiff as in the 
optical trap experiments in [13]). 

An elastic force F(x) = k(x — xq) represents the sim- 
plest and most generic cc-dependent force. Whereas for a 
confinement of fixed length or a constant force, the MT 
length x was the only stochastic variable, the force F(x) 
itself is now coupled to x and becomes stochastic as well. 



Therefore, not only are the MT length distributions of in- 
terest but also the maximal and average polymerization 
forces which are generated during MT growth. 



A. Vanishing rescue rate 

We first discuss growth in the absence of rescue events, 
uj r = 0. This situation corresponds to optical trap exper- 
iments [11, 13], which are performed on short time scales 
and no rescue events are observed. In a state of growth 
the MT grows against the elastic obstacle with velocity 
Vj r [f{x)\ and f(x) increases. For simplicity we suppress 
the x-dependency in the notation in the following. At a 
maximal polymerization force / max , the MT undergoes 
a catastrophe and starts to shrink back to zero and the 
dynamics stop due to missing rescue events. No steady 
state is reached. Since switching to the state of shrinkage 
is a stochastic process, the maximal polymerization force 
/max is a stochastic quantity which fluctuates around its 
average value. We calculate the average maximal poly- 
merization force (/max) within a mean field approach. 
Here (...) denotes an ensemble average over many real- 
izations of the growth experiment. 

Because no steady state is reached in the absence of 
rescue events, we have to use a dynamical mean field 
approach, which is based on the fact that the MT growth 
velocity dx/dt = f+(/) is related to the time evolution of 
the force by df/dt = (k/F )dx/dt. In mean field theory, 
this results in the following equation of motion for (/), 



id) = !« + «/», 



(29) 



where we used the mean field approximation (v + (j)) « 
v + ((/))■ With the initial condition (f)(0) = we find a 
time evolution 



(/)(*) 



In 



(l-a; on /w off )e * /r +w on /w off (30) 



/stall + In [1 - exp(-i/r)] 



with 
10 2 ... 



a characteristic time scale r 
10 4 s for k w 10- 5 ...10- 7 N/m. 



(31) 

= F Q /dkuj oS w 
For long times 
i > r, cq. (30) approaches the dimensionless stall force 
(/) = /stall, see eq. (14), which is the maximal polymer- 
ization force in the absence of catastrophes. The approx- 
imation (31) holds for Won/^off' 3> 1- 

MT growth is ended, however, by a catastrophe, and 
the average time spent in the growing state is t = 
l/w c ((/ max )). Together with eq. (30), this gives a self- 
consistent mean field equation for the maximal polymer- 
ization force (/max), 

(/max) = In [(1 - Won/Woff) e^/^/—^ + . 

(32) 

The maximal polymerization force (/ max ) is always 
smaller than the stall force /stall as can be seen from 
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FIG. 5: Average maximal polymerization force (/max) for 
an elastic obstacle and in the absence of rescues as a function 
of Won for different values of k = 10" 5 N/m(A), 10" 6 N/m(B) 
and 1CP 7 N/m(»). Data points represent results from simula- 
tions, solid lines are solutions of eq. (32). Error bars repre- 
sent the standard deviation of the stochastic quantity (/max). 
Dashed line: dimensionless stall force /stall = In (uj on /u! s). 



B. Non-zero rescue rate 

For a non-zero rescue rate uj r , phases of growth, in 
which f(x) increases and which last l/w c (/) on average, 
are ended by catastrophes which arc followed by phases 
of shrinkage. Shrinking phases last l/u) r on average, and 
during shrinkage the elastic obstacle is relaxed and f(x) 
decreases. After rescue, the MT switches back to a state 
of growth. In contrast to the case without rescue events, 
the system can attain a steady state. In this steady state, 
the average length loss during shrinkage, V-/ui r , equals 
the average length gain during growth, v+(f)/uj c {f) 1 and 
the MT oscillates around a time-averaged stall length (x) , 
which is directly related to the time-averaged polymer- 
ization force by (/) = (k/Fo)((x) — xq)- In the following, 
the steady state dynamics and the average polymeriza- 
tion force are characterized. We start with an analysis 
of the full master equations focusing on the stationary 
state followed by a dynamical mean field theory, which 
can also be applied to dilution experiments. 

In the presence of a x-dependent force f(x), the master 
equations for the time evolution of p+^—(x,t) become 



eqs. (30,31). Since u> on /u> ff ^ cj c t ^> 1 for realistic force 
and parameter values, eq. (32) can be approximated by 



111 



Woff™ c ((/ max )) 

/stall - In [™ c ((/ max )) 



(33) 



For a catastrophe rate increasing exponentially above the 
characteristic force Fq, eq. (15), we find 



i 



ax) ~ „ In ( ,„s 

2 V F o uj c {0) 



-,dk 



(34) 



i.e., the maximal polymerization force grows logarithmi- 
cally in uj on (note that the catastrophe rate in the absence 
of force decreases as w c (0) cx l/u; on [26]), see Fig. 5 for 
k = 10~ 5 N/m. Within a slightly different catastrophe 
model obtained from experimental data and discussed in 
section VII, this logarithmic dependence can be shown 
exactly. 

Fig. 5 shows (/max) as a function of u> on . Analyti- 
cal results from eq. (32) agree with numerical findings 
from stochastic simulations. The maximal polymeriza- 
tion force (/max) increases with increasing k, see eq. (34), 
but it remains smaller than the stall force /stall- Stochas- 
tic simulations show considerable fluctuations of /max, 
which arc caused by broad and exponentially decaying 
probability distributions for / ma x and which we quantify 
by measuring the standard deviation (/ max ) — (/max) • 
For increasing k, probability distributions become more 
narrow and mean field results approach the simulation 

results for (/max)- 



d t p+(x,t) = —U) c (x)p+(x,t) +U r P-(x,t) 

- d x (v+(x)p + (x,t)) (35) 
d t p-{x,t) = u c (x)p + {x,t) - w r p-(x,t) + V-d x p-(x,t), 

(36) 

which differ from eqs. (1) and (2) by the x-dependence of 
growth velocity and catastrophe rate. Both growth ve- 
locity v + (x) = v + [f(x)] and catastrophe rate oj c (x) = 
ui c {v+[f(x)]} become x-dependent via their force- 
dependence. Therefore, also the force-dependent length 
parameter A(/) from eq. (26) becomes x-dependent via 
its force-dependence, X(x) = A[/(x)]. Eqs. (1) and 
(2) are supplemented by reflecting boundary conditions 
v + (0)p + (0,t) = v_p_(0,i) at x = 0, similar to eq. (3). 

For the steady state, eqs. (35) and (36) are solved on 
the half-space x > with reflecting boundary conditions 
at x = 0, and we can calculate the overall MT length 
distribution P{x) = p+{x) +p_(x) explicitly, 



P(x) = N 1 + 



v + (x) 
with a normalization 



e *oM(o) exp 



dx'/X(x') 



N- 1 = 



dx[l + — n - 

V v +l x ) 



(37) 

o x /\(0) e f* dx'/\(x') ^ ^ 



where X(x) = A(/ = 0) in the force-free region x < xq and 
X(x) = X[f(x)] for x > xq and, likewise, v+(x) = v + (f = 
0) for x < xq and v+(x) = v+[f(x)] for x > Xq. This 
implies e Xo ^ x< -°^e^o dx ^ x( - x ' — e x / x (°) an d, thus, a simple 
exponential dependence of P(x) for x < Xq. A similar 
OPDF has been found for dynamic MTs in the presence 
of MT end-tracking molecular motors [36] . 
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With increasing length x, also the force f(x) increases 
and, thus, v + [f(x)] decreases and to c [f{x)] grows expo- 
nentially. If x becomes sufficiently large that the con- 
dition A -1 [/(a;)] < holds, the distribution P(x) starts 
to decrease exponentially. In this length regime the MT 
undergoes a catastrophe with high probability. Because 
the distribution always decreases exponentially for suffi- 
ciently large x, a single MT growing against an clastic 
obstacle is always in the regime of bounded growth re- 
gardless of how large the values of w on and u r are chosen. 
This behavior is a result of the linearly increasing force, 
which gives rise to arbitrarily large forces for increasing x 
in contrast to growth under constant or zero force, where 
a MT can either be in a phase of bounded or unbounded 
growth as mentioned above. 

The behavior is also in contrast to length distribu- 
tions in confinement between fixed rigid walls, where we 
found a transition between exponentially decreasing and 
increasing length distributions: The elastic obstacle typ- 
ically leads to a non-monotonic length distribution with 
a maximum in the region x > xq. (as long as the on- 
ratc uj on and rescue rate w r are sufficiently large and 
the obstacle stiffness k sufficiently small). While res- 
cue events (and an exponential decrease in the growth 
velocity v + [f(x)]) cause P{x) to increase exponentially 
for small MT length, catastrophes are responsible for 
an exponential decrease for large x. The interplay be- 
tween rescues and catastrophes gives rise to strongly lo- 
calized probability distributions with a maximum. Figs. 
6 (a-d) show the steady state distribution P(x) obtained 
from eq. (37) for different values of Lu on and to r . We 
chose k = 10~ 7 N/m and xq = 10 /im. In the steady 
state, a stable length distribution with a well defined 
average length (x) = J Q °° P{x)xdx is maintained al- 
though the MT is still subject to dynamic instability. 
The length distributions drop to zero for large x, where 
A _1 (.t)~— lu c (x)/v+(x) and uj c (x)/v + (x) increases expo- 
nentially with increasing force. 

The most probable MT length x rnp maximizes the sta- 
tionary length distribution (37). Because u_ 3> v+(x) 
and using the approximation of an exponentially decreas- 
ing growth velocity, v + \f(x)} ps v + (0)e~^ x \ which is 
valid for ui on ^> w ff (see eq. (13)), we obtain a condition 
A _1 (a; TO p) = -d x f(x mp ) = -k/F or 

V+(fmp)Ur - V-U)c(fmp) = -(k/F )V-V + (f mp ) (39) 

for the corresponding most probable force f mp = 
(k/F )(x mp - x ). 

For an exponentially increasing catastrophe rate above 
the characteristic force F , eq. (15), we find 



/, 



V-LO C (0) 



1 



F u r 



(40) 



We can distinguish two limits: (i) For a soft obstacle with 
kv-/Foui r <C 1 the most probable force f mp is identical 
to the critical force f c for MT dynamics under constant 
force, see (28), because the right hand side in the condi- 
tion (39) for f mp can be neglected and we exactly recover 



condition (27) for f c . The most probable MT length thus 
"self-organizes" into a "critical" state with f mp ps f c , 
and a MT pushing against a soft clastic obstacle gen- 
erates the same force as if growing against a constant 
force. This force grows logarithmically in the on-rate 
o; on and the rescue rate u r . (ii) For a stiff obstacle with 
kv-/FoLd r 1, on the other hand, the most probable 
force is larger than the critical force, f mp 3> / c , and the 
MT growing against a stiff obstacle generates a higher 
force. This limit can also be realized for vanishing rescue 
rate u) r , and for kv—/FoU) r > 1 wc indeed recover the 
maximal pushing force in the absence of rescue events, 
i-e. f mp * (/max) from eq. (34) with v + (0) ps w on d. This 
force grows logarithmically in the on-rate u> OIL . Further- 
more, if f mp becomes negative for small on-rates and res- 
cue rates (leading to A _1 (0) < —k/Fo, see eq. (40)) the 
stationary length distribution has no maximum, see for 
example Figs. 6(a,b) at the lowest on-rates. 

With respect to the MT's ability to generate force the 
two limits can be interpreted also in the following way: 
Fq is the characteristic force above which the catastro- 
phe rate increases exponentially. For kv-/Fou> r <C 1, the 
average length loss during a period of shrinkage, V—/uj r , 
is much smaller than the length Fo/k, which is the dis- 
placement x — xo of the elastic obstacle under the charac- 
teristic force Fq. Therefore, the MT tip always remains 
in the region x > xq under the influence of the force for 
a soft obstacle with kV-/FaOJ r -C 1, whereas it typically 
shrinks back into the force-free region x < xq before the 
next rescue event for a stiff obstacle kv-/FoOj r ^> 1. The 
force generation by the MT can only be enhanced by res- 
cue events if rescue takes place under force in the regime 
x > xq. Therefore, we find an increased polymerization 
force f mp ps f c > (/max) as compared to the force / max 
without rescue events discussed in the previous section 
only in the limit kv-/Fou> r <C 1, i.e., for a soft obstacle or 
sufficiently large rescue rate. In the limit kv-/Fouj r 3> 1 
of a stiff obstacle, the MT only generates the same force 
as in the absence of rescues, f mp ps (/max)- 

By comparing the condition (27) or v+(f c ) = 
v~.uj c (f c )/uj r for the critical force f c , the condition (39) or 

V+{f m p) = V-U c (f mp )/U) r {l + fcl'_/F ) < V-U) c (f mp )/(J r 

for the most probable force f mp , and the condition 
u+(/staii) = for the stall force, see eq. (14), it follows 
that 



fc — fmp ^ / 



/stall 



(41) 



i.e., force generated against an elastic obstacle is be- 
tween critical and stall force but typically well below the 
stall force, which is the maximal polymerization force in 
the absence of catastrophes. Therefore, the stall length 
Zstali = (F /k) In (w on /w off ) + xq is always much larger 
than the most probable MT length x mp at the maximum 
of the stationary length distribution, see Fig. 6(a). This 
shows that the dynamic instability reduces the typical 
MT length significantly compared to simple polymeriza- 
tion kinetics. 

In order to quantify the width of the stationary dis- 
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tribution P(x) we expand the exponential in (37) up to 



second order about the maximum at x r 
first expand A -1 (a;) up to first order: 



To do so we 



^0 



v + (x mp )oj T + n_w c (x mp ) 



(■>■ 



V + (x mp )V- 

(42) 

where we used v + [f(x)] ~ v + (0)e f( x \ which is valid 
for w on 3> w g- (see eq. (13)), and where we approxi- 
mated the catastrophe rate by an exponential uj c [f{x)\ ss 
w c (0)e^^^ according to eq. (15) resulting in Lu' c [f(x)] ~ 
ku c [f(x)]/ F$. The prime denotes a derivative with re- 
spect to the length x. Using the expansion (42) in eq. 
(37), we obtain an approximately Gaussian length distri- 
bution 



P(x) « N 1 + 



v+(x) 



e x /\(0) x 



exp 



with a width 



a 2 = 



(■f, 



Xq 



2g 2 



exp 



(x 



2<T 2 



(43) 



k 



V+(x mp )uJ r + V-UJ C 
\ 2 / 1 



i 



(44) 



where we used the saddle point condition (39) in the 
last approximation and the exponential approximations 
v+[f(x)] « w+COje-^W and w c [/(s)J « w c (0)e^). 
Again we have to distinguish the two limits of soft and 
stiff obstacles: (i) For a soft obstacle with kv-/FoU) r <C 1 
we find a 2 rj FoV~/2ku) r . This shows that the width of 
the length distribution decreases with increasing w r but 
is roughly independent of the on-rate w on , as can also 
be seen in the series of simulation results shown in Figs. 
6. Closer inspection of the simulation results shows that 
the width of the stationary length distribution P{x) is 
slightly decreasing with the on-rate ui on . (ii) For a stiff 
obstacle with ki)—/ Fquj t 3> 1, on the other hand, we find 
a 2 w (F /k) 2 , which only depends on obstacle stiffness. 
All in all, a 2 is monotonously decreasing for increasing 
stiffness k. 

For a soft obstacle kv-/Foui r <C 1, high rescue rates 
thus lead to a sharply peaked length distribution P(x) 
and suppress fluctuations of the MT length around x = 



and we expect (x) 



to a very good approxi- 



mation. This property of a sharp maximum in P{x) will 
make the mean field approximation that is discussed in 
the next section very accurate. 

If the obstacle stiffness k is increased the most probable 
MT length x mp = xq + f mp Fo/k approaches Xo, and a 
considerable probability weight is shifted to MT lengths x 
below xq (see Fig. 6 (e)). The average length approaches 
and finally drops below xq. This signals that the force 
generated by the MT is no longer sufficient to push the 
obstacle out of its equilibrium position xq- The obstacle 



now serves as a fixed rigid boundary and P(x) approaches 
the results eq. (21) and (22). The dynamics of a single 
MT within confinement can therefore be seen as a special 
case of the dynamics in the presence of an elastic obstacle, 
i.e., for small w on and w r or for large spring constants k. 

So far we have quantified the generated force by the 
most probable force f mp - The generated force can also 
be quantified by the average steady-state force (/) = 
Jn°° f(x)P(x)dx. Using the stationary distribution (37) 
with normalization (38) we can calculate (/); results are 
shown in Fig. 7 in comparison with the most probable 
force f m p, which is determined numerically from the max- 
imum of P{x), and the stall force / s t a ii m the absence of 
dynamic instability from eq. (14). For (/), there is ex- 
cellent agreement with stochastic simulations over the 
complete range of parameter values. The results clearly 
show that the dynamic instability reduces the ability to 
generate polymerization forces since, even for large val- 
ues of w on and w r , the average force (/) is always smaller 
than the stall force. Nevertheless forces up to F ~ 1.5 Fo 
can be obtained in the steady state for realistic param- 
eter values. Comparing (/) and f mp we find (/) < f mp , 
and both forces become identical, (/) rj fmp, in the limit 
of large rescue rates or a soft obstacle kv-/F uj r <C 1, 
where also the length distributions P{x) become sharply 
peaked, sec Fig. 6. Comparing different combinations of 
w on and u) r and the corresponding forces, one finds that 
the influence of the on-rate w on on force generation is 
more significant than the influence of the rescue rate u> r . 
For w on = 100 s -1 , a four fold increase of the rescue rate 
w r gives rise to an increase of (/) by a factor of ~ 1.5, 
while for u> r = 0.1s -1 , a four fold increase of the on-rate 
w on results in an amplification of the force (/) by a fac- 
tor of ~ 9. These results can be explained within a mean 
field theory presented in the next section. 



C. Mean field approach (non-zero rescue rate) 

In the following, we show that we can reproduce many 
of the results for the average polymerization force (/) 
for non-zero rescue rate using a simplified mean field 
approach. Using the mean field approach, we can also 
address the time evolution of the average force (/), for 
example, in dilution experiments. Since the switching 
between the two states of growth is a stochastic pro- 
cess, the length x and the force f(x) are stochastic vari- 
ables. Therefore, the velocity of growth v+[f(x)] and the 
catastrophe rate w+[/(a;)] also become stochastic vari- 
ables which, in the steady state, fluctuate around their 
average values. Within the mean field approach we ne- 
glect these fluctuations and use (v + [f(x)]) = v + ((f)) 
and (w+[f(x)]) = w+ ((/)). In the mean field approxi- 
mation, the average time in the growing state is given 
by l/w c ((/)) and the average growth velocity is v + ((/)). 
The average time in a shrinking state is l/w r . There- 
fore, the mean field probabilities to find the MT grow- 
ing or shrinking are p+ = w r /[w r + w c ((/))] and p_ = 



11 



, 0.1 2 

g 0.08 
^0.04 



(a) w r =0.03s _1 
a; s taii(w on = 25s _1 ) 




25 50 75 100 




(c) 




25 s -1 5 


0s- 1 75S-H00S- 1 



25 50 75 100 125 



25 50 75 100 125 




100 150 



0.5 
0.4 




K = J.U IN/ m I c 1 


0.3 




W on = SOS" 1 

ui t = 0.05 s- 1 


0.2 
0.1 




j 


fc = 10- 6 N/m 
A fc = 10_- 7 N/m 



40 60 x[fj,m] 



FIG. 6: Stationary MT length distribution P(x) in the steady state for growth against an elastic obstacle with uj oi1 — 
25s" 1 , 50s" 1 , 75s" 1 , 100s" 1 and different values of uj r . We set k = 10" 7 N/m and x = 10" 5 m. (a) oj r = 0.03s" 1 . (b) 
uj r = 0.05 s" 1 . (c) uj t = 0.1 s" 1 . (d) uj r = 0.2 s" 1 . Dashed line represents xq. In picture (a) the stall length x B for uj on = 25 s" 1 , 
obtained from simple polymerization kinetics, is indicated by an arrow, (e): P(x) for uj oi1 = 50s" 1 , uj t = 0.05s" 1 and different 
values of the spring constant k. 
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FIG. 7: Average steady state force {/} as a func- 
tion of uj on for growth against an elastic obstacle with 
tor = O.O3s" 1 (a),O.O5s" 1 (B),O.ls~ 1 (0),O.2s" 1 («) and k = 
10" 7 N/m. Solid lines: (/) = J °° f(x)P(x)dx with P(x) given 
by eqs. (37) and (38). Dashed lines: (/) calculated from mean 
field equation (46). Dotted lines: most probable force f mp , 
measured in simulations, for ui r = 0.03 s" 1 and ui r = 0.05 s" 1 . 
Also shown is the dimensionless stall force /stall obtained from 
simple polymerization kinetics (14). 



uj c ((f))/[u} r + oj c ((f))], respectively. This results in the 
following mean field average velocity v of a single MT 
under force: 



v((f)) 



UJ r + W c ((/)) 



(45) 



In the steady state the barrier is pushed so far that (/) 
stalls the MT. We require v((f)) = and obtain the 
condition 



v+({f))u r = u_w c «/» 



(46) 



for the stationary state. This condition corresponds to 
a force, where the average length gain during growth, 
v +((f))/ UJ c ((f)), equals the average length loss during 
shrinking, From the mean field equation (46), 

the average steady state force, (/) can be calculated as 
a function of u r and w on - The average length (x) can 
be obtained from the relation (/) = (k/Fo)((x) — xq). 
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FIG. 8: Average steady state force (/) as a func- 
tion of Won for growth against an elastic obstacle with 
uj r = O.O3s" 1 (Q),O.O5s" 1 (B),O.ls" 1 (0),O.2s" 1 («) and k = 
10" 5 N/m. Solid lines: (/) = J °° f(x)P{x)dx with P(x) 
given by eqs. (37) and (38). Dashed lines from bottom 
to top: (/) calculated from mean field equation (46) for 
uj r = 0.03 s" 1 , 0.05 s" 1 , 0.1 s" 1 and 0.2s" 1 . 



Results obtained from the mean field equation (46) match 
numerical results from stochastic simulations very well as 
shown in Fig. 7. 

The mean field condition (46) is identical to the condi- 
tion (27) for the critical force f c for MT dynamics under 
constant force such that 



(/> = fc 



(47) 



which can be interpreted as "self-organization" of the 
average MT length or the average force to the "critical" 
state. Therefore, the curves presented in Fig. 7 for (/) 
are identical to the curves shown in Fig. 4 (b) for f c . 

This also allows us to take over the results we derived 
for the critical constant force f c . Using the approxi- 
mation of an exponentially decreasing growth velocity, 
v + [f(x)] w v+(0)e~^ x \ which is valid for uj on ^ w Q ff 
(see eq. (13)), and an exponentially increasing catastro- 
phe rate above the characteristic force Fq, eq. (15), we 
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find 



2 \v-uj c (0) 



(48) 



which is identical to the result (28) for f c . 

Comparing with the stall force and the most probable 
force, we use relation (41) and find 
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(/> = fc < fmp < /stall- (49) 

In the limit of a soft obstacle, fcv_/i*oCJ r -C 1, the average 
force (/) approaches the most probable force (/) ps f mp , 
whereas the mean field average force (/) is always smaller 
than the stall force /stall in the absence of dynamic in- 
stability from eq. (14). 

Finally, we discuss the limits of validity of the mean 
field approximation. The mean field approximation is 
based on the existence of a pronounced maximum in 
the stationary MT length distribution P(x), which con- 
tains most of the weight of the probability density P(x). 
It breaks down if this maximum broadens or vanishes, 
such that a considerable amount of probability density 
is shifted below xq into the regime of force-free growth. 
Then the MT typically shrinks into the force-free region 
x < xq during phases of shrinkage such that the grow- 
ing phase explores the whole range of forces starting from 
/ = up to / > (/), and the approximation of a constant 
average force / ps (/) during growth is no longer fulfilled. 
For small spring constants k or large values of u) r , the 
length distribution P(x) assumes a Gaussian shape with 
width it, see eqs. (43) and (44). When k is increased for 
a fixed combination of iu on and u r , the average length (x) 
approaches xq as (x) — xq oc 1/fc, whereas the width a of 
the length distribution only decreases as a oc 1/ yfk in the 
regime of a soft obstacle kv-/FQiu r <C 1, as can be seen 
from eq. (44). Therefore, an increasing amount of proba- 
bility density is shifted below x , where no force is acting 
on the MT ensemble (see Figs. 6) (a) and (e)). The mean 
field approximation is only valid for spring constants k 
which fulfill (x) —xo S> er/2 for given parameters w on and 
io r . With (/) = (k/Fo)((x) — xq) this is equivalent to a 
condition 

ka If 2F uj r \~ 1/2 , . 

according to eq. (44). This condition can only be fulfilled 
in the limit of a soft obstacle with kv-/FoU) r <C 1. For 
the validity of the mean field approximation we therefore 
recover the condition that the average length loss during 
a period of shrinkage, V-/ui r , is much smaller than the 
typical displacement F /k of the elastic obstacle under 
the characteristic force F . Then the MT tip always re- 
mains in the region x > x under the influence of the 
force. 



FIG. 9: (a): Average force (/)(<) as a function of time for 
k — 10~ 7 N/m, uj r — 0.05 s _1 and different values of u> on . 
Symbols: time dependent average force (f)(t) measured in 
simulations. Solid lines: time dependent average force trajec- 
tory calculated from eq. (51). (b): Average force {f)(t) as a 
function of time for k — 10~ 7 N/m, uj on — 50 s -1 and different 
values of uj r . Symbols: time dependent average force (f)(t) 
measured in simulations. Solid lines: time dependent average 
force trajectory calculated from eq. (51). 

D. Dynamics and dilution experiments 

Within the mean field approach we can also derive an 
analytical time evolution of the average time-dependent 
force (f)(t). The time evolution is based on eq. (45), 
which gives a mean field approximation for the average 
MT velocity v((f)) as a function of the average force. 
On the other hand, the average MT growth velocity is 
related to the time derivative of the average force by 

> " JE5 W " To° m (51) 

Using eq. (45) for v((f)), this gives a mean field equation 
of motion for (/) similar to eq. (29) in the absence of 
rescue events. Integrating this equation numerically we 
obtain mean field trajectories for the average force (/) (t) 
as a function of time t. Figs. 9 shows such trajectories 
for k = 10~ 7 N/m and a initial condition (f)(0) = at 
t = 0. Also shown in Figs. 9 are results from stochas- 
tic simulations, which show excellent agreement with the 
mean field trajectories. 

We now address the question of how fast a single MT 
responds to external changes of one of its growth pa- 
rameters. Here we focus on fast dilution of the tubulin 
concentration, which is directly related to the tubulin on- 
rate w on . In vivo tubulin concentration can be changed 
by tubulin binding proteins like stathmin [37] , while in in 
vitro experiments, the tubulin concentration can be di- 
luted within seconds [38] . In the following we give a mean 
field estimate of the typical time scale, which governs the 
return dynamics of the MT back to a new steady state 
after the tubulin on-rate is suddenly decreased. In the 
initial steady state the average velocity v((f)i) vanishes 
and the average polymerization force (/).; (and, thus, the 
average length (x)i) can be calculated from the condition 
v + ((f)i)uj r = v-ui c ((f)i), cf. eq. (46), for a given com- 
bination of w on and uj r . If uj on is suddenly decreased 
this leads to a sudden decrease in the growth velocity to 
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FIG. 10: Average force (f)(t) as a function of time t. Symbols 
are results obtained from simulations. We set k = 10~ 7 N/m, 
Lo r = 0.05s -1 and w n = 75s . At t = 20000s, u> on is diluted 
down to ujon = 50 s -1 . Solid line represents a fit with an ex- 
ponential decay (53) to the simulated data with fit parameter 
Td — 1762 s. Dashed lines indicate the average force in the 
initial state /, before dilution and in the new final state // 
after dilution. 



v+(f) < «+(/) and an increase of the catastrophe rate 
to to c (f) > u c (f), resulting in a negative average veloc- 
ity ««/)) = [«+«/> H - u-w c «/»]/[a> P +w c ((/))] < 
according to eq. (45). Consequently, the MT starts to 
shrink with an average velocity < 0. This relaxes 

the force from the elastic obstacle, i.e., (f)(t) starts to 
decrease from the initial value /, = (/),. With decreasing 
average force (/) (t), the average growth velocity v((f) (t)) 
increases again (because v+ increases and uj c decreases) 
until the steady state condition v+((f) f)to r = V-ui c ((f) f) 
holds again and a new steady state force (/)/ < (/); is 
reached (s. Fig. 10). 

The relaxation dynamics to the new steady state af- 
ter tubulin dilution is therefore governed by the average 
velocity «((/)) given by eq. (45). To extract a charac- 
teristic relaxation time scale, we expand the average ve- 
locity v((f)) to first order around the final steady-state 
polymerization force // = (/)/, which is the solution of 
eq. (46) with uj r and the decreased tubulin on-rate w on , 
which takes its dilution value. Using = one finds 

in first order 



««/» 



V + (ff)w r + V-U' c {ff) 



UJ r + Ul c (ff) 



((f) -ff) (52) 



where the prime denotes the derivative with respect to 
the force. In the last approximation we used the mean 
field condition eq. (46) and v + [f(x)] ps v + (0)e~^ x \ 
which is valid for uj on ^> cj g (see eq. (13)). This ex- 
pansion is only valid for average forces close to the new 
average polymerization force ff. Using this expansion, 
the time evolution (51) of the average force after dilution 
exhibits an exponential decay 



(f)(t) = ff + (h-ff)e- t/Td 
with a characteristic dilution time scale 

F u r +oj c (ff) F uj r + uj c (ff) 



Td 



v +(ff)Ur + V-U' c (ff) k 2V-U) c (ff) 



(53) 



(54) 



where we approximated the catastrophe rate by an expo- 
nential u) c [f(x)] m ui c (0)ef( x ) according to eq. (15), and 



we used the mean field condition eq. (46). In the limit 
w c(//) 3> £*> r , i.e., at forces ff ^> 1, we obtain the sim- 
ple result Td ss Fo/2v-k. In general, the relaxation time 
Td is proportional to the square a 2 of the width of the 
stationary distribution, cf. eq. (44): A narrow length dis- 
tribution gives rise to fast relaxation to the new average 
force. 



VII. EXPERIMENTAL CATASTROPHE MODEL 

So far we have employed the catastrophe rate derived 
by Flyvbjerg et ai, to which we will refer as w C; Fiyv m 
the following. This expression for the catastrophe rate 
was based on theoretical calculations of the inverse pas- 
sage time to a state with a vanishing GTP-cap, see eq. 
(10). In order to investigate the robustness of our results 
with respect to changes of the catastrophe model, we now 
investigate an alternative expression for the catastrophe 
rate that has been obtained from experimental results. 
Throughout this section, we focus on the third confine- 
ment scenario of an elastic obstacle, and we compare re- 
sults from the two different catastrophe models for zero 
rescue rate u r = and non-zero rescue rate ui r > 0. In 
addition, we restrict the comparison to mean field results, 
since numerical and stochastic calculations match mean 
field results well over the complete range of parameters 
(see Sec. VI). 

Experimentally, it has been found that the average 
time (t+) spent in a growing state is a linear function of 
the growth velocity v + [39]. The force-dependent catas- 
trophe rate is then given by 



^c,Jans(/) 



av+(f) + b 



(55) 



and 



with constant coefficients a = 1.38 
b = 20 s. At «+(/) = 0, w c , Jans (/) = 0.05 s- 1 and 
for v + (f) = —b/a, the catastrophe rate w c ,jans(/) di- 
verges. This is in contrast to the theoretical model, where 
w c ,Fiyv(/) is finite for all v+{f). Also w C: j ans (/) increases 
exponentially for forces F > Fq or / > 1. This com- 
mon feature is essential and lead to similar results for 
both catastrophe models. In Fig. 11, both catastrophe 
rates are shown as a function of the dimensionless force 
/. The catastrophe model (55) is based on experimental 
data and, thus, is phenomcnological. It assumes neither 
a purely chemical model, as in the model by Flyvbjerg et 
al. , nor a chemo- mechanical model in the sense of "struc- 
tural plasticity" [24] . 



A. Vanishing rescue rate 

We start with the case uj t = without rescue events, 
and we calculate the average maximal polymerization 
force within the experimental catastrophe model using 
the self-consistent mean field eq. (32), which holds inde- 
pendently of the choice of catastrophe model (see Sec. 
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FIG. 11: Catastrophe rate ui c (f) as a function of force / for 
ojon = 50 s~ and oj on = 75 s _1 . Solid lines: uj Cy F\ yv from the 
catastrophe model by Flyvbjerg et al. Dashed lines: w c ,Jans 
from the experimental catastrophe model by Janson et al. 
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FIG. 12: (a): Average maximal polymerization force (/max) 
as a function of ui on and w r = for k = 10~ 7 N/m, 10~ 6 N/m 
and 10 _7 N/m (top to bottom) . Dotted line: dimensionless 
stall force /stall- Solid lines: (/max) obtained from tJ Cl Fiy V (eq. 
(32)). Dashed lines: (/max) obtained from w c ,jans (eq. (56)). 
(b): Average steady state force (/) as a function of td on . k = 
10~ 7 N/m, bj r = 0.2 s _1 (top) and ui r = 0.03 s _1 (bottom). 
Solid lines: (/) obtained from a; Ci pi yv (eq. (46)). Dashed lines: 
(/) obtained from o; Ci j ana (eq. (60)). Dotted line: dimension- 
less stall force /stall- 



B. Non-zero rescue rate 



VIA). As for the catastrophe by Flyvbjerg et al., we 
have w c .j ans r 3> 1 for realistic parameter values and 
v+((f)) < —b/a, and eq. (32) can be solved explicitly 
for (/ m ax) in this limit. We find an average maximal 
polymerization force 



(fn 



In 



(A 2 



B) 1/2 -A 



(56) 



with 



(cJon/^off - l)adw ff - (<Jon/Wf - 1)6 - T 



2t 



B = 



(w on /u oS - l)adu 



Since cu on /u! ff ^> 1, eq. (56) can be approximated by 

(/max) ~ hi (w on /w max ) (57) 



with 



0J,, 



2tu! q h 



2 1 1/2 

(adu) ff — b) + 4adcj ffT — [adw Q ff — 6] 

(58) 

For realistic parameter values, we have r 3> adw ff > 6, 
and recover the expression (34) derived using the Flyvb- 
jerg catastrophe model: 



(fn 



1 



hi 



1 



hi 



-.dk 



F cj, 



c.Jans 



(0) 



(59) 



In Fig. 12 (a), (/max) as obtained from eq. (32) with the 
Flyvbjerg catastrophe model and eq. (56) with the ex- 
perimental catastrophe model are shown as a function of 
cj on . Results match qualitatively and quantitatively well, 
although they are obtained from two different catastro- 
phe models. The maximal polymerization force (/ max ) 
always remains smaller than the stall force /stall- 



Now we compare both catastrophe models for a non- 
zero rescue rate, and we calculate the average steady 
state force. For the experimental catastrophe rate (55), 
the mean field equation (46) can be solved explicitly, and 
the average steady-state force (/) is given by 



(/) = ln(a; on /cj av ) , 



(60) 



ith 



(2 \ 1 / 2 

(i) ~id +UJ ° s (61) 

Again (/) < / staU since w av > w off . Fig. 12 (b) show 
(/) as a function of w on . For realistic parameter values, 
we have V-/u> r 3> b 2 /a and (iL/wrarf 2 ) 1 ' 2 3> u> Q s, and 
recover the expression (48) derived using the Flyvbjerg 
catastrophe model: 



</>-m 



1 / v+(0)oj r 

2 V^-^c.JansCO) 



(62) 



In Fig. 12 (b), results for (/) from both catastrophe 
models are shown as a function of on-rate w on . The 
average steady state force obtained from w Ci Fiyv is al- 
ways slightly larger than (/) obtained from cj C) j ans , since 
<^e,jans(/) > <^e,Fiyv(/) for forces smaller than or compa- 
rable to F . Otherwise, both results agree qualitatively 
and quantitatively well. 



VIII. FORCE- VELOCITY RELATION 

Finally, we discuss the influence of the force-velocity 
relation on the MT dynamics. We restrict our analysis to 
mean field results obtained for the third scenario, i.e., the 
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clastic obstacle. A change in the force-velocity relation 
directly modifies the velocity of growth «+(/), but it also 
affects the catastrophe rate uj c (v + (f)), which are both 
crucial parts of the MT dynamics. In the following, we 
employ a more general form of the force- velocity relation, 
which is consistent with thermodynamic constraints, and 
we show that our results are robust with respect to this 
generalization. 

In their investigation of experimental data Kolomeisky 
et al. used a generalized growth velocity 

v+{f,0) = d{w on exp(-0/) - w off exp[(l - 8)f}}, (63) 

which depends on a dimensionlcss "load distribution fac- 
tor" 9 [40]. The load distribution factor 9 £ [0, 1] deter- 
mines whether the on- or off-rates are affected by external 
force, while keeping the ratio of overall on- and off-rate 
unaffected. Under force both the tubulin on-ratc uj on 
and the tubulin off-rate w ff now acquire an additional 
Boltzmann-likc factor. For 9 — 1, we obtain again v + (f) 
as given by eq. (13). The dimensionlcss stall force is un- 
affected by 9 and is still given by /stall = m ( w on/ w off)- 



A. Vanishing rescue rate 

We use the generalized force- velocity relation v+(f, 9) 
given by eq. (63) and the catastrophe rate w c ,Fiyv(/) m 
order to calculate the average maximal polymerization 
force (/max) from the self-consistent mean field eq. (32). 
In Fig. 13 (a), (/ max ) is shown as a function of the load 
distribution factor 9 for k = 10~ 5 N/m and different val- 
ues of Won- At 9 = 1, the maximal force (/max) equals the 
maximal polymerization force obtained with v+(f) from 
eq. (13). With decreasing 9, (/max) increases but remains 
below the dimensionlcss stall force. The growth velocity 
v + (/, 9) increases with decreasing 9 for a fixed force / 
and, therefore, the maximal polymerization force (/max) 
increases. For high tubulin on-rates uj on = 75 — 100 s -1 
and small 9 rj 0, . . . , 0.2, the maximal polymerization 
force (/max) approaches the dimensionlcss stall force. 



B. Non-zero rescue rate 

For non-zero rescue rate, the average steady state force 
(/) is calculated from the mean field eq. (46), where we 
use the force- velocity relation v+(f,0) (eq. 63) and the 
catastrophe rate cj c ,Fiyv(/)- in Fig. 13 (b), results for 
(/) are shown as a function of 9 for k = 10 _7 N/m, 
u> r = 0.05 s" 1 and different values of uj on . At 9 = 1, 
(/) equals the average steady state force obtained with a 
velocity v+(f) taken from eq. (13). The average steady 
state force (/) increases with decreasing 9, as explained 
above. For high tubulin on-rates uj on = 75 — 100 s" 1 and 
small 9 « 0, . . . , 0.2 , also the average steady state force 
(/) again approaches the dimensionlcss stall force but 
remains smaller. 



(a) (b) 






FIG. 13: (a): Solid lines: Average maximal polymerization 
force (/max) as a function of 9 for k = 10~ 5 N/m and different 
values of iu on . Dashed line: Dimensionless stall force /stall for 
ojon = 100s -1 . (b): Solid lines: Average steady state force 
(/) as a function of 9 for k — 10~ 7 N/m, uj r = 0.05 s _1 and 
different values of u OIX . Dashed line: Dimensionless stall force 
/stall for w on = 100 s -1 . 



IX. DISCUSSION AND CONCLUSION 

We studied MT dynamics in three different confining 
scenarios: (i) confinement by fixed rigid walls, (ii) an 
open system under constant force, and (iii) MT growth 
against an elastic obstacle with a force that depends lin- 
early on MT length. These three scenarios represent 
generic confinement scenarios in living cells or geome- 
tries, which can be realized experimentally in vitro. For 
all three scenarios, we are able to quantify the MT length 
distributions. In scenario (iii) of an elastic obstacle, 
stochastic MT growth also gives rise to a stochastic force. 
For this model, we also quantify the average polymeriza- 
tion force generated by the MT in the presence of the 
dynamic instability. 

The parameter A, see (4) and (26), governs the MT 
length distributions in confinement by fixed rigid walls, 
and under a constant force. For confinement by rigid 
walls we introduced a realistic model for wall-induced 
catastrophes. There is a transition from exponentially in- 
creasing to exponentially decreasing length distributions 
if A changes sign. The average MT length is increasing for 
increasing on-rate and increasing rescue rate, as shown 
in Figs. 2. Wall-induced catastrophes lead to an overall 
increase in the average catastrophe frequency, which we 
quantify within the model. 

For MT growth under a constant force, there exists a 
transition between bounded and unbounded growth as in 
the absence of force. This transition takes place where 
the parameter A(/) changes sign. Under force, the tran- 
sition to unbounded growth is shifted to higher on-rates 
or higher rescue rates and determines a critical force / c , 
see Figs. 4. 

MT growth under a MT length-dependent linear elastic 
force allows for regulation of the generated polymeriza- 
tion force by experimentally accessible parameters such 
as the on-rate or the rescue rate. The force is no longer 
fixed but a stochastically fluctuating quantity because 
the MT length is a stochastic quantity. For zero rescue 
rate, i.e., in the absence of rescue events, we find that 
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the average maximal polymerization force (/max) before a 
catastrophe depends logarithmically on the tubulin con- 
centration and is always smaller than the stall force in 
the absence of dynamic instability as shown in Fig. 5. 

For a non-zero rescue rate, we find a steady state length 
distribution, which becomes increasingly sharply peaked 
for increasing rescue rate and is tightly controlled by 
microtubule growth parameters, see Figs. 6. Interest- 
ingly, the average microtubule length self-organizes such 
that the average steady state polymerization force (/) 
equals the critical force for the boundary of bounded and 
unbounded growth, (/) = f c . Because of the sharply 
peaked MT length distribution, the average polymeriza- 
tion force (/) can be calculated rather accurately within 
a mean field approach as can be seen in Figs. 7 and 8. 
The average polymerization force is always smaller than 
the stall force in the absence of dynamic instability. 

Within this mean field approach, we can also describe 
the dynamics of the average force, see Figs. 9. This might 
be useful in modeling dilution experiments, where the re- 
sponse to sudden changes in the on-rate is probed. For 



this type of experiment, we estimate typical polymeriza- 
tion force relaxation times. 

Finally, we show that our findings are robust against 
changes of the catastrophe model (Figs. 12) as long as the 
catastrophe rate increases exponentially above a char- 
acteristic force and that results are also robust against 
variations of the relation between force and polymeriza- 
tion velocity in the growing phase (Figs. 13), which are 
obtained by introducing a load distribution factor. 
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